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Abstract 

In this paper, we derived Lorentz covariant quantum Liouville equation for 
the density operator which describes the relativistic quantum information pro- 
cessing from Tomonaga-Schwinger equation and an exact formal solution for 
the reduced-density-operator is obtained using the projector operator tech- 
nique and the functional calculus. When all the members of the family of the 
hypersurfaces become flat hyperplanes, it is shown that our results agree with 
those of non-relativistic case which is valid only in some specified reference 
frame. To show that our new formulation can be applied to practical prob- 
lems, we derived the polarization of the vacuum in quantum electrodynamics 
up to the second order. The formulation presented in this work is general 
and could be applied to related fields such as quantum electrodynamics and 
relativistic statistical mechanics. 
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Recently, there has been growing interest in the relativistic formulation J7| of quantum 
operations for possible near future applications to relativistic quantum information process- 
ing such as teleportation ||, entanglement-enhanced communication ||, and quantum clock 
synchronization |K|, [ [□]] . 

In the non-relativistic case, the key element for studying quantum information processing 
is the density operator of a quantum register which is derived from the solution of a quantum 
Liouville equation (QLE) fll2| , [ |13|1 for the total system including an environment. The QLE 
is an integro-differential equation and it is in general nontrivial to obtain the solution of the 
form 

P^p' = i\p], (i) 

where p is the reduced density operator of the quantum register and £ is the superoperator 
describing the evolution of p by the quantum information processing. In the previous works, 
we have employed a time-convolutionless reduced-density-operator formalism to model quan- 
tum devices [frj| and noisy quantum channels JIB] . 



The first step toward the relativistic quantum information theory would be the formu- 
lation of Lorentz covariant QLE and the derivation of the reduced-density-operator which 
is a solution of the covariant QLE. The goal of this paper is to derive Lorentz covariant 
quantum Liouville equation which describes the relativistic quantum information processing 
and obtain a formal solution for the reduced-density-operator pertaining to the system (or 
electrons) part alone. 

It is well known that neither the non-relativistic Schrodinger equation nor the QLE is 
Lorentz covariant. As a result, it is expected that the usual non-relativistic definition of 
the reduced-density-operator and its functionals such as quantum entropy have no invariant 
meaning in special relativity. Another conceptual barrier for the relativistic treatment of 
quantum information processing is the difference of the role played by the wave fields and the 
state vectors in the quantum field theory. In non-relativistic quantum mechanics both the 
wave function and the state vector in Hilbert space give the probability amplitude which 



can be used to define conserved positive probability densities or density matrix. On the 
other hands, in relativistic quantum field theory, covariant wave fields are not probability 
amplitude at all, but operators which create or destroy particles in spanned by states defined 



as containing definite numbers of particles or antiparticles in each normal mode |T7|]. The 
role of the fields is to make the interaction or S-matrix satisfy the Lorentz invariance and 
the cluster decomposition principle. The information of the particle states is contained in 
the state vectors of the Hilbert space spanned by states containing 0, 1, 2, • • • particles as in 
the case of non-relativistic quantum mechanics. So it seems like that one needs to obtain 
the covariant equation of motion for the state vector and derive the covariant QLE out of 
it. 



Some time ago, Tomonaga |18| and Schwinger derived a covariant equation of motion 
for the quantum state vector in terms of the functional derivative, known as Tomonaga- 
Schwinger (T-S) equation, 



i- 



H int (x)V[a], (2) 



5a(x) 

in the interaction picture. Here x is a space-time four-vector, a is the spacelike hypersurface, 
ty[cr] is the state vector which is a functional of a, Hi nt (x) = l~LintYPa{ x )\ is the interaction 
Hamiltonian density which is a functional of quantum field y9 Q [x], and is the Lorentz 
invariant functional derivative [ 20| . The functional derivative of \I/[<t] is defined as 



—LA = i im U 3 

Sa(x) duo 

where 5u is an infinitesimal four- dimensional volume between two hypersurfaces a and a'. 
The formal solution of equation (0) is given by 

tt[a]=W[<7,<7oM<7b], (4) 

where the generalized transformational functional satisfies the T-S equation 

i 5 -^^=H int (x)U[a,a ] (5) 



with the boundary condition W[<7o,o"o] = 1. The generalized transformation functional 



U[a, uq] is a unitary operator. We also have ]R 



W [<x, cr J U [a, cr J, 



(6) 



5cr(a;) ' 5<r(x) 

from the unitary condition. Throughout the paper, we assume ft = c = 1. The expectation 
value of some field variable _F(x) becomes 

< F(x) > = (V[a],F(x)V[o]) 

= trace(F(x)^[a]^[a]) 

= trace(F(x)pr[cr]). (7) 
^From equation (0), we notice that the total density operator pr[(r] can be written as f21 



22 



p r [cr] = ^[aj^^cr] 

= U[<J,a ]9 [(Tq] & [a ]U- 1 [a, a ] . 



(8) 



Then, 



5cr(x) 



£{Wfo<7o]*to]*Vo]W~Wo]} 

00" 



<5W[o-,cto] 



5a(x 
-i[H int (x),p T [a]] 

-iC(x)p T [cr), 



a,a \,p T [a\ 



(9) 



where C(x) is the Liouville superoperator. Since equation @ describes the Lorentz covariant 
equation of motion for the total density operator, we denote it as the covariant quantum 
Liouville equation (CQLE). Note that the Liouville superoperator is not an operator in the 
Hilbert space of state vectors but a linear operator in the Hilbert-Schmidt space of density 
matrices |T7J. Here pr\?\ contains the information for the total system, for example, an 
interacting spin-- massive particles and photons in the case of quantum electrodynamics 
(QED). 



In order to extract the information of the system or the electrons alone, it is convenient 



to use the projection operators [|12|], P3] , f24}| that decompose the total system by eliminating 
the degrees of freedom for the environment, say, the photon field in the case of QED. The 
information of the system is then contained in the reduced-density-operator p[a] which is 
defined as 

p[a] = tr B p T [a] 

= tr B Vp T [cr], (10) 

where the projection operator V and Q are defined as VX = p^tr^X), Q = 1 — V, for any 
covariant dynamical variable X, ps is the density matrix for the quantum environment at a~o 
and trs indicates a partial trace over the quantum environment. The projection operators 
satisfy the operator identities V 2 = V, Q 2 = Q, VQ = QV = and [j^,V] = [i, Q] = 
0. Furthermore, we would like to note that (^^j) -1 = J d 4 x pOU , and the system and the 



environment are decoupled at <jq- We also note that the projection operators V and Q 
are functionals of the initial hypersurface a (^ a for all x) and unless otherwise specified, 
we will omit the functional argument. However, one needs to keep track of the functional 
argument especially in the four- dimensional integration. 

The CQLE (|9p can be decomposed into two coupled equation for Ppr[c] and Qpt[o]: 
5 



5a(x 
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-Vp T [a] = -iVC{x)Vp T [a\ - tVC(x)Qp T [a] } (11a) 
■ Qp T [a] = -iQC(x)Qp T [a} - iQ£(x)Vp T H (lib) 



5a(x) 

In order to obtain the formal solution, we solve first eq. ( |llb[ ) using the integrating 
factor. Let h[a] be an integrating factor such that 

h[a}{-~-Qp T [cT]+iQC(x)QpT[(j]} = -ih[a)QC(x)Qp[a] 
da{x) 

' {h[a]Q P [a]}. (12) 



Sa(x 

Then, ^jfy = ih[a]QC(x)Q and we obtain 
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h[a] =T c exp{t f d*x'Q£{x')Q}. (13) 

/From eq. (|TJ), 

J (TO 

= -i ^ d 4 xH[a(x),a(x , )}QC(x')Vp T [(T(x% (14) 

J (TO 

where we assume that Pt[&] is decoupled when a = cr and 

H[(j{x),a{x')]=Texp{-i ' d A x"QC{x")Q}. (15) 



Here T and T c are time-ordering and anti-time ordering operators, respectively, and 
H[a(x), c{x')] is the projected propagator. In order to derive the convolutionless equation 
of motion, we define the retarded propagator G R [a(x), a (x')] such that 



G R [a{x),a{x')]=T c exp{i F * d 4 x"C(x")} (16) 

which satisfies 

p T [<r ] = G R [a, (T ] p[u] ■ (17) 

Then, 

Qp[a] = -i r d 4 x'H{a(x),a(x')]QC(x')VG R {a(x),a(x')]pT{a(x)] 



<T() 



d 4 x'H[a{x),a(x')}Q£(x')VG R [a(x),a(x')}Vp T W(x)} 



(T() 



i f 7 d 4 x'H[a(x),a(x')]Q£(x')VG R [a(x),a(x')]Qp T [a(x)} (18) 

J (TO 



and 



Qp T [a] = {9[a} - l}Vp T H (19) 



where 



0- 1 [ a ] = l + t f d 4 x'H[a(x),a(x')}Q£.(x , )VG R [a(x),a(x% (20) 

J (TO 



Once the solution for Qpt[o~] is obtained, it is substituted for the equation for Vpt[o]- Then, 
after, some mathematical manipulations, we obtain using the integrating factor technique 
again, 

Vp T [a] = W- l [a, a }U s [a, a ]Vp T [cr }, (21) 



or 



p[a] = tr B {W- l [a, a ]U s [a, a ]p B }p[a ], (22) 



where 



W[a, (7 ] = 1 + i r d 4 x'U s [a(x),(i(x')]VC(x'){9[a(x')] - l}VG R [a{x), a(x')]9[a(x)], (23) 

J (T 

and 

U s [a,a ] = Texp{-i f d 4 x'VC(x')V}. (24) 

J (JO 

Here W s [cr, cr ] is the generalized transformation functional or the propagator for the reduced 
system. 

It is remarkable to note when hypersurfaces ctq and all the members of the family {a} 



are hyperplane flat surfaces parametrized by t = constant [20|, then the transformation 



functional such as U s [a(x),a(x')} can be written as U s (t,t'). Then, if we set to = 0, 

W(t,0) = W{t) 

= 1 + / d 4 x'U s {t,t')V£{x',t){6(t') - l}VG R (t,t')9(t) 
Jo 

= 1 + / dsU s (t, s)tr B [C(s){9{s) - l}p B )tr B [G R (t, s)6(t)p B ], (25) 
Jo 

with C(s) = J d 3 x'(x', s). As a result, the covariant forms of equations ( pi]) to ( p4|) become 
reduced to those of the non-relativistic case which is valid only in some specified reference 
frame given by equations (18) to (25) of reference JT5|j . 

By comparing, equations (|I|) and (|2"2"D, the covariant superoperator for the relativistic 
quantum operation S[a, o"o] can be written as 



£[*, (7 ] = trgiW-^a, a ]l( s [a, cx ]p B }. (26) 



So far all our results are exact and the equations ([H]) to ( pq) would be the key steps in 
the analysis of relativistic quantum information processing. Apart from describing quantum 
information processing, QLE and reduced-density-operator have been essential in solving 
various quantum optics and non-Markovian optical problems in the non-relativistic domain 
||12|| - ||14)| . So it might be interesting to extend this approach to revisit relativistic quan- 
tum electrodynamics problems, which were solved relying on renormalization procedures in 
field theory, using the covariant form of quantum Liouville equation. On the other hand, 
relativistic thermodynamics or statistical mechanics look like an area where the knowledge 
of the density operator or the reduced-density-operator might come in handy provided the 
ambiguity of the temperature concept in special relativity is resolved. We believe our formal- 
ism is general and could be applied to related fields such as QED and relativistic statistical 
mechanics. As a matter of fact, these related fields would also play an important role in 
relativistic quantum information processing because these processes would cause the deco- 
herence as in the non-relativistic case. 

To show how to apply the formalism we developed to practical problems, we give a 
derivation of the polarization of the vacuum by an external field starting from the equation 
([22|). The Hamiltonian for the coupling between an electron and electromagnetic fields is 
given by 

H(x) = -j^A^x), (27) 

where j^(x) and A^x) are current and electromagnetic 4- vector potential operators, re- 
spectively. Then from eq. (^2|) , the reduced density operator up to the first order in C 
becomes 

p0)[ a ] = tr B {l-i r d A x f PC(x')V}p[a ] 

= (1-i f d 4 x'tr B (C(x')p B ))p[a ]. (28) 

J err, 
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If we set the initial hypersurface be the flat surface &q = — oo, p[— oo] = p and A^(x) 
tr B (A fJi (x)p B ) which is a classical external field, we get 



pW[<r] = p + i r d 4 x%{x')A,(x'),p ). (29) 



oo 



The polarization of the vacuum is the expectation value of j^(x), computed for the state of 
the system as modefied by the external electromagnetic field ||25|| , [P6fl and is given by 



= trQ^po) + i [ d 4 x'tr{\j v (x')A u (x'),pQ\]^x)} 

J — oo 

= i [ d A x'tr{[] ll (x),i v (x)}po}A u (x') 



a 1 

H) 



15p -nV,(x) + --- (30) 



where £r(- • •) is the trace over the electron states and (• • -}o is the expectation value for the 
electron fields. Here Jp(x) is the external current generating the electromagnetic field, fco = 
moc/h, a = e 2 /4irhc, D 2 = <9 M <9 M , and tuq is the electron mass [^5] | 26[| . Eq. (|30|) describes the 



vacuum polarization due to the external electromagnetic fields in quantum electrodynamics. 
We proceed to derive the second order correction to the vacuum polarization (j^x))^ 2 ' . The 
second order correction term to the reduced density operator Ap^[a] becomes 



ApW[a] = -i f d A x'tr B {(W- l [a,-oo]l( s [a,-oo]Y 2) p B }po 

J -co 

= -2 r [X) d A x' r {X) d A x"tr B {C{x')C{x")p B )p 

J —oo J —oo 

= -2 r {x) d 4 x' r {x \ 4 x"{-(A^x>)A„(x")) j,(x')Ux")p 

J —oo J —oo 

+ (A^x')A u (x")) j u (x")p j^x') + (A u (x")A^x')) j^x') P oUx") 
-(A v (x")A^x')) p Ux")Ux')}- (31) 



Then 



Q,{x))V =trQ,{x)ApV[a]) 



-2 r (x) d 4 x f r [x,) dV{-(4(x')i,(x"))o(^(x),i(x')]i(x"))o 

J —oo J —oo 

+(A u (x")A IJ ,(x')) Q„ (x") [j x (x) , %{x') ] ) 0} 

-2 f W dV r (X ' ] rfV([i,(x'), i,(x"))o([jA(x), j M (x')]i(x"))o 



—00 •' —00 



-2i f {X) d A x' r (X ' ) dVV^(^'-^)(bA(a;),j M (x')]j,(x // ))o 
J —00 J— 00 

-2i r W dV f ^ d 4 x"([j,(x),i(^)]i(^))o^(^' - ar"), (32) 



00 j —00 



where T>(x) is the invariant function defined by eq.(2.17) of reference Above result can 
be further simplified by using that 2H 



d A x"V(x' — x ,, )jf i (x n ) 

00 

00 

d 4 x"e(x',x")V(x' -x")j^x") 

-00 

00 _ 

d 4 x"V(x' - x'^j^x") 

-00 

= -25A^) (33) 

Here SA^x) is the 4- vector potential induced by the polarization of the vacuum or the 
reaction of the virtual electron-positron coupling. Then the vacuum polarization up to the 
second order interaction becomes 

Qp{x)) = 1 r iX) d 4 x'^^x),Ux')]) Q {A v {x') + A8A v {x')} (34) 

J —00 

The knowledge of vacuum polarization would be important in understanding the decoher- 
ence process in the relativistic domain. At this stage, we would like to leave the detailed 
calculations of the second and higher order corrections for the future work. 

In summary, we have derived Lorentz covariant quantum Liouville equations for the den- 
sity operator in functional of hypersurface from T-S equation and obtained formal solution 
for the reduced-density-operator which is also in covariant form using the projection op- 
erator technique and the functional calculus. When all the members of the family of the 
hypersurfaces become flat hyperplanes, our results agree with those of the non-relativistic 
case. We have shown that our new formalism can be applied to the practical cases such as 

10 



the vacuum polarization. Our formulation is exact and general so it could be applied not 
only to the relativistic quantum information processing but also to the related fields such 
as QED, field theory and relativistic statistical mechanics. 
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